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We propose a technique to prepare coherent superpositions of two nondegenerate quantum states 
in a three-state ladder system, driven by two simultaneous fields near resonance with an intermediate 
state. The technique, of potential application to enhancement of nonlinear processes, uses adiabatic 
r*j , passage assisted by dynamic Stark shifts induced by a third laser field. The method offers significant 

advantages over alternative techniques: (i) it does not require laser pulses of specific shape and 
duration and (ii) it requires less intense fields than schemes based on two-photon excitation with non- 
^ ' resonant intermediate states. We discuss possible experimental implementation for enhancement of 

00 ' 

CNj ' frequency conversion in mercury atoms. 
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O PACS numbers: 42.50.Hz, 32.80.-t 
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' I. INTRODUCTION 

Efficient mechanisms to generate coherent superpositions of quantum states are central to a rich variety of applica- 



X 



tions in modern quantum physics. Quantum logic gates, i.e. the key components of a quantum computer [1], rely on 
superpositions of two degenerate quantum states as qubits [2, 3]. Quantum tunnelling and localization in a double well 
potential can be controlled [4] by techniques that require preparation of coherent superpositions. Numerous authors 
have noted that nonlinear optical processes, e.g. resonantly enhanced frequency mixing in atomic vapours, can be 
significantly improved [5-10] by preparing the nonlinear optical medium in a coherent superposition of nondegener- 
ate quantum states. Whereas numerous techniques exist to prepare coherent superpositions of degenerate states, as 
needed for qubits in quantum computing, nonlinear optics still hold challenges. For example, frequency conversion to 
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short wavelength radiation involves high-lying states, requiring excitation by multi-photon transitions. 

Procedures based on adiabatic passage driven by coherent interactions [11] provide reliable and robust tools for the 
creation of superpositions. Unlike diabatic techniques, which rely on precise control of individual pulses, adiabatic 
processes are insensitive to small variations of pulse duration and peak intensity. High-intensity laser systems, com- 
monly used for efficient frequency conversion, exhibit fluctuations in intensity and other parameters. Thus robust 
techniques are most appropriate to support nonlinear optical processes driven by the high-intensity lasers. 

In what follows we will discuss adiabatic passage in a three-state system involving two near-resonant laser pulses. 
We assume that initially the atom is in the ground state 1 of energy E\ and that the first laser field, of frequency ui\, 
links this with an excited state 2 of energy E 2} while the second field, of frequency u>2, links this state with a final 
target state 3 of energy E 3 . We assume that the two fields, though near resonant with the specified transitions, are far 
from resonance with any other transition. The relative ordering of the third-state energy E 3 is not significant; we shall 
assume that it lies above E 2 , as is appropriate for application to nonlinear optics, so that the linkages form a ladder. 
Figure 1 shows this linkage pattern, along with spectroscopic labels appropriate to implementation in mercury. 

4 

n 3 . 



(7 1 So) 



3- 

A 5 



—J. A, 



ft; 



ftz 



(6 3 Pi) 2 



fti 



(6%) 1 



FIG. 1: Three-state system in ladder linkage, interacting with two pulsed laser fields, indicating the detunings Aj, the Rabi 
frequencies Q.j and the Stark shift As. Dotted lines show transitions for four-wave mixing, through a nonresonant state 4. To 
the left are the spectroscopic labels appropriate to implementation with mercury. 



Our modeling is based upon the three-state time-dependent Schrddingcr equation in the rotating wave approxima- 
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tion (RWA) [12], for which the Hamiltonian is 



ifli(t) 







H(t) = h 



ifii(t) A 2 



*fi 2 (t) 



(1) 



in 2 (t) A 2 + A 3 -A s (i) 



The elements of this RWA Hamiltonian are the detunings 



A 2 = (E 2 -Ei)/h-uj-L 



A 3 = (£3 - E 2 )/h-oj 2 . 



(2) 



and the pulsed Rabi frequencies fii(i) = d\ 2 E\{i)l% and ri 2 (i) = d<2.%£i{i)lh that parametrize the laser-induced 
excitation interaction with the dipole-transition moments d nm of the transition between states n and m and the 
electric field envelopes £i(i) and £ 2 (t). A third laser, nonresonant with either of these transitions, induces dynamic 
Stark shifts of the energies. For the ladder-like linkage the shift is expected to be largest for state 3; we denote that 
shift as Ag and we neglect shifts of the other states. 
In the RWA the statevector has the expansion 



where ip n is a unit vector representing quantum state n. Our principal objective is to create, for times t later than 
the conclusion of a pulse sequence at tj, an equal-probability superposition of states tpi and ip^, specifically 



where ip is a time independent phase defining the relative sign of the superposition. This particular superposition, 
with equal probabilities P n (t) = \C n {t)\ 2 = \{Tp n Y$>{t))\ 2 of the two constituent nondegenerate states, provides the 
basis of a technique referred to as "nonlinear optics at maximum coherence" [9], that substantially improves the 
efficiency of nonlinear frequency conversion processes. The following section reviews how the efficiency of four- wave 
mixing is enhanced by preparing the atomic medium in the coherent superposition of eqn. (4) with equal probability 
amplitudes. We review in section III two techniques recently suggested to generate coherent superpositions of this 
form: (i) fractionally-completed stimulated Raman adiabatic passage (F-STIRAP) and (ii) half-completed Stark- 



*(t) - Ci(t)Vi + CtWe-^ik + C 3 (i)e- j(wi+aJ2)t V<3 



(3) 



t>t f 



(4) 
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chirped rapid adiabatic passage (half-SCRAP). Section IV describes our proposed method, Stark-assisted coherent 
superposition (SACS), in which two near-resonant fields, accompanied by a laser-induced AC Stark shift, produce 
adiabatic passage into the desired superposition. Section V illustrates different techniques by presenting simulations 
of excitation of mercury vapor, a medium of significant interest for applications in frequency conversion. 

II. ENHANCEMENT OF FOUR- WAVE MIXING BY MAXIMUM COHERENCE 

The starting point for nonlinear optical phenomena is the wave equation describing the propagation of the electric 
field E through matter. The effects of matter are incorporated into a polarization P that serves as an inhomogeneous 
term in the wave equation, 

E(i,r)= Mo ^P(i,r). (5) 

We idealize the medium as a uniform distribution of identical motionless atoms, of number density Af. Then the 
polarization is the Af time the expectation value (d(i,r)) of the single-atom dipole moment: 

P(i,r)=JV<d(t,r)>. (6) 

In the following we idealize the laser fields as plane waves traveling along the z axis. We suppose that there are several 
such fields, each characterized by a common polarization vector e but different carrier frequency u) k - We express the 
resulting electric vector as 

E(i,r) = ^2e£ k (t')cos(uj k t' + (j> k ) (7) 
k 

where t' =t — z/c. The presence of the laser fields causes changes of the material. With these changes come changes in 
the individual dipole expectation values and hence of the polarization P. We will discuss a four- wave mixing process 
in an atomic gas, as are typical for the generation of short-wavelength radiation, i.e. vacuum- or extreme-ultraviolet. 
For application to four-wave mixing in mercury, the three states depicted in Fig. 1 arc: (1) 6 1 So, (2) 6 3 Piand (3) 
7%. 

For the four- wave mixing processes, two fields at frequencies u)\ and L02 create an induced dipole moment. This 
dipole moment, together with a third field of frequency W3 far from resonance with any state V>4, combine to produce 
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a dipole moment that varies at the frequency u 4 = u\ + UJ2 + ^3- Following the definition of the expectation value of 
the dipole moment (d(t')) = (^(t')\d\^/(t')) , it can be expressed as a function of the probability amplitude 



(d(O)-e = 2Re [c? (f)C 2 tf)d Vi e iuit ' + C;{t')C 3 (t')d 23 e W2t ' 



(8) 



The last term serves as the source of an electric field at the frequency 0J4. For a detuning A4 = (£4 — E 3 )/h — ui 3 
much larger than the Rabi frequency O3, one can adiabatically eliminate the state ip4 [12]. The amplitude C^t') is 
then given by 



2A4 



2A4 



(9) 



When f^4 is small compared to £l 3 , as it is initially, the source term of the electric field at frequency UJ4 is 



r 4 (t') = C* 3 (m(t')£ 3 (t') d ^e^' 



(10) 



where Vj(t), for j=(l,2,3,4), is defined from the expression 



P(i') = 2eRc 



(11) 



These successive equations show that the induced polarization, and hence the efficiency of any subsequent nonlinear 
optical process, depends on the product C3 (t')Ci(t'), i.e. on the coherence, established between the ground and 
the excited state. Laser fields typically exhibit inhomogeneities in time, space or phase. If there is a dynamical 
contribution to the phase, these inhomogeneities will, when averaged, tend to reduce any coherence to nearly zero. 
To generate the desired radiation at the frequency uj 4 it is thus crucial that the relative phase of the components 
Ci(t) and C 3 (t) in the superposition state of eqn. (4) does not depend on the time integrated Rabi frequencies, as 
dynamical phases do. If this condition is fulfilled, the source term of eqn. (10) reaches a maximum value when the 
atomic coherence is maximum, i.e. when |Ci(t)| = | C3 (t) | . The polarization, the efficiency of four- wave mixing and 
thus the intensity of the generated radiation at frequency 104 will be enhanced, if the atomic medium is prepared in 
maximum coherence (4). 
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III. ADIABATIC PREPARATION OF COHERENT SUPERPOSITIONS 

As will be discussed below, to produce this superposition we make use of adiabatic states and adiabatic 

eigenvalues, Afc(t), defined as instantaneous solutions to the eigenvalue equation of the RWA Hamiltonian 



*T(t)*fc(t) - ftA fc (t)$fc(i). 



(12) 



Prior to describing the proposed technique, we comment on two common adiabatic techniques, first intended for 
complete population transfer but later applied to the preparation of superposition states. Figure 2 illustrates the 
three schemes. 
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FIG. 2: Three-level system in ladder-type configuration for the specific adiabatic excitation processes: (a) STIRAP, (b) F- 
STIRAP (c) half-SCRAP, (d) SACS. Below the linkage patterns are indicators of the time dependence of the pulses. 



A. STIRAP and F-STIRAP 

Stimulated Raman adiabatic passage (STIRAP) [13, 14] is well established as an appropriate tool to drive complete 
population transfer in a three-state lambda-type system [11]. With suitable modification [15-20], the technique 
provides a means of creating superpositions of initial and final states. Essentially the evolution is adiabatically 
stopped just when half of the population is transferred. We refer to this as fractional STIRAP (F-STIRAP). 

The basic STIRAP procedure is based on a lambda linkage in which E 3 < E 2l i.e. a stimulated Raman process (see 
Fig. 2a). It works also in a ladder linkage if the lifetime of the uppermost state is much longer than the interaction 



time, i.e. the pulse duration - which we will assume in the following. The process relies on interaction of the 
system with two laser pulses, termed Stokes and pump, which produce overlapping but not coincident interactions, 
parameterized by Rabi frequencies, Cls{t) and Qp(t). The Stokes pulse, linking states 2 and 3, precedes the pump 
pulse, linking states 1 and 2. Adiabatic evolution requires that the time integral of the Rabi frequencies is much 
larger than 1. Under these conditions (and two-photon resonance, E% — E\ = hu>\ + hu>2) the statevector remains at 
all times aligned with one particular adiabatic state, the dark state or population trapping state, 



*o(t) 



n s (t) 
o 

-e-^n P {t) 



(13) 



where is the difference between Stokes and pump phases. At the beginning of the interaction in F-STIRAP, i.e. 
when the Stokes is strong and the pump pulse is negligible, the adiabatic state <&o(i) coincides with the initial state 
V>i- STIRAP proceeds to completion when the pump field is strong and the Stokes field is negligible, i.e. at the end 
of the interaction. The adiabatic state <I>o(t) is then aligned with ^3, and complete population transfer has occurred. 

By contrast, in F-STIRAP the pulses need to be terminated simultaneously (see Fig. 2b). Once the pump pulse 
has produced the desired fractional transfer, both pulses must diminish simultaneously while maintaining a fixed 
ratio of the two Rabi frequencies, £lp(t)/tts(t) — tana. As with STIRAP, it is necessary to maintain the two-photon 
resonance condition A2 + A3 = 0. When the mixing angle a remains constant as the pulses diminish, the result is 
the statevector 



*( Q )(t) = cos a Vi - sine 



,-i(3 -i(wp+v s )t 



'3, 



t>t f . 



(14) 



(Note the phases that appear here, a consequence of the steady rotation of the coordinates of the adiabatic states.) 
The desired superposition of the ground and target states, eqn. (4), requires a = 7r/4, i.e. equal Rabi frequencies as 
the pulse sequence concludes. Because F-STIRAP, in contrast to regular STIRAP, demands laser pulses of specific 
temporal shape and duration, it is difficult to implement experimentally in a system having a ladder-type configuration. 
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B. 



Half-SCRAP 



Two-photon excitations in ladder systems, involving absorption of two photons from one radiation field, permit 
population transfer to high-lying excited states without the need for lasers with very short wavelength. However, 
because such two-photon interactions involve nonresonant couplings to intermediate states (which can be summarized 
in a virtual state), strong two-photon excitation is inevitably accompanied by dynamic Stark shifts. Such shifts, when 
produced by nonresonant fields, can be put to good use as a means of sweeping the two-photon detuning through 
resonance and thereby inducing Stark chirped rapid adiabatic passage (SCRAP) [21, 22]. Like STIRAP, SCRAP was 
initially developed as a technique for complete population transfer. Subsequently, it has been used to create transient 
superpositions for use in enhancing frequency conversion [10]. A modification of SCRAP has been suggested for the 
preparation of persistant superposition states, the so-called half-SCRAP process [23]. 

The half-SCRAP technique (see Fig. 2c) employs an intense pump laser field that couples, e.g. via a two-photon 
transition, a ground state if>i and a target state tp 3 . There is no resonance with an intermediate state. The interaction 
that links states ipi and ip 3 is represented by an effective two-photon Rabi frequency Q e ff(t) proportional to the 
intensity of the pump laser. The pump field induces dynamic Stark shifts of the energies E n . An additional Stark- 
shifting laser pulse induces further energy shifts, which should be larger than the Stark shifts induced by the pump 
laser. Experimentally this is easy to implement, required a strong, fixed-frequency Stark-shifting laser is available. 
Typically the polarizability of the upper state greatly exceeds the polarizability of the ground state, thus only the 
Stark shifts of the energy E 3 need be considered. The resulting detuning of the pump laser from two-photon resonance 
is 



where lo p is the pump-field carrier frequency, HAp(t) is the (small) energy shift induced by the pump field and ftAg(t) 
the (large) shift induced by the Stark-shifting field. Under appropriate conditions the combination of the two pulses 
will produce, apart from an overall phase factor, the superposition 



Aeff(t) = {Es -E x )jh- 2uj v - A s (t) - A P (t) 



(15) 



* e (t) = cos 6>Vi - sinfle-^e- 2 ^* 



V>3 



t>t f 



(16) 



where the mixing angle 9, defined by the equation tan2#(i) = Q e s(t) / A e s(t) , is here evaluated at the time tf, and 
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where (3 is the phase of the pump field. 

The half-SCRAP process requires that the following conditions be met: (i) the Stark-shifting pulse must precede 
the pump pulse, (ii) the adiabatic condition, \8(t)\ <C A e s(t) 2 + Q e g(t) 2 , must be fulfilled and (iii) at the end of the 
interaction the effective two-photon Rabi frequency f2 e ff(i/) must be much larger than the effective detuning A e g(tf). 

Although the half-SCRAP technique relies on an adiabatic process, and hence offers the robustness of any adiabatic 
process, it has two practical disadvantages. As noted, the pump-induced two-photon coupling f2 e ff(i) is accompanied 
by a dynamical Stark shift included in A ff (t) . These two interactions are both proportional to the intensity of the 
pump pulse. The pump laser-induced dynamic Stark shift has to be compensated by an additional static two-photon 
detuning, such that at the end of the interaction fl c ff(tf) ^> \A e s(tf)\. In contrast to techniques utilizing one-photon 
transitions, half-SCRAP as described above requires a pump laser of large intensity to drive a two-photon transition 
sufficiently strongly. 



IV. STARK-ASSISTED COHERENT SUPERPOSITION (SACS) 

A. Basic principles 

We propose an alternative technique to prepare the coherent superposition of eqn. (4). The mechanism, termed 
Stark-assisted coherent superposition (SACS), is based on adiabatic passage and uses the dark state of eqn. (13). In 
contrast to STIRAP and F-STIRAP two simultaneous, rather than delayed radiation fields with a similar temporal 
shape are applied. As in SCRAP, a non-resonant Stark-shifting pulse varies the energy E 3 by —hAs(t) (see Fig. 2d). 

Our starting point is the RWA Hamiltonian of eqn. (1). In the following we assume the two Rabi frequencies to 
be equal, Qi(t) — 2 (i) = Q(t), and we require two-photon resonance, A 2 + A 3 = 0. Without loss of generality we 
assume that the Stark shift is negative (Ag(i) > 0). Thus we deal with the RWA Hamiltonian 



H(t) = h 



JrQ(f) 

ifi(t) A 2 ift(t) 
*fi(i) -A s {t) 



(17) 



Because we focus on adiabatic evolution, the statevector ^(t) for this system should at all times be aligned with an 
adiabatic state <&&(£) of this Hamiltonian. The eigenvalues Afc(£) and eigenvectors $&(£) are readily obtained for any 
time t by numerical means. In two important limiting cases they have simple properties that underly our proposed 
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method. When the interaction is absent, i.e. SI = 0, but a Stark shift (and positive A 2 ) is present, the three adiabatic 
states align with the diabatic states ip n . The eigenvalues and eigenstates are then 



A_=-A s , $_(t) 



A = 0, $ (t) 



A+ = A 2 , $+(t) 





1 

1 





1 





1p3, 



^1 



(18) 



(19) 



(20) 



When, in addition, no Stark shift is present, then the adiabatic states $_ and <j> are degenerate and can be taken as 
any superposition of states ipi and ipz] we take these to be the choices presented here. 

In the other limiting case, when no Stark shift is present, A5 = 0, the interaction O moves the adiabatic states 
away from alignment with the diabatic states. The eigenvalues are then 



A_ 
A 
A+ 



1 
2 

0, 

1 
2 



a 2 - Va 2 + 2n 2 



A 2 + VA 2 + 2fi 2 



(21) 
(22) 
(23) 



Of particular interest is the null-eigenvalue adiabatic state; i.e. the "dark state" superposition 



7! [ ^ ~ H 



(24) 



The proposed excitation process can be understood by observing the three adiabatic eigenvalues A& (t) as they vary 
due to changes of the RWA Hamiltonian. Figure 3 displays these three energies, evaluated numerically, as a function 
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of the two parameters Q and As that define this Hamiltonian. Similar figures have earlier been used to explain 
adiabatic processes; the topology of the surfaces can exhibit conical intersections and avoided crossings [24-26]. 




FIG. 3: (Color online) Surfaces of eigenenergies as function of the Rabi frequency Q and of the Stark shift As (in units of 
A2). The path represented by the (blue) curve on the intermediate surface gives values of the couple (Q, As) that connect the 
initial state ip\ to the superposition of eqn. (4). 

The two-photon resonance condition ensures that initially, when no interaction is present, two of the adiabatic 
states, $0 and 3>- share the eigenvalue zero. A large dot on the figure indicates this point in parameter space (the 
initial system point that characterizes the system in the parameter space as the initial statevector characterizes it in 
the time space). At this time the remaining eigenstate, <£>+, aligned with the bare state ifc, has the eigenvalue A2, 
also denoted by a dot in the figure. When O = the three adiabatic states align uniquely with the diabatic basis 
states, tp n . 

The SACS process can be described by a path on a surface in the parametric space of Fig. 3. At the beginning no 
fields are present; the statevector is aligned with state tpi and with the adiabatic state $0. Subsequently the detuning 
As grows, while the Rabi frequency remains null. In the vertical plane at f2 = the energy line with constant value 
is associated with the adiabatic state $o ; while the constant value A2 belongs to adiabatic state $+. The adiabatic 
state $_ has an energy that varies linearly with the Stark shift As- Because there is no interaction to alter the 
composition of the adiabatic states, the statevector remains aligned with the initial state ip\ during the growth of the 
Stark-shifting pulse. Once As has reached a satisfactory value, the interaction fields O are applied. Thereafter the 
system point moves, on the energy surface of $01 away from the plane at Q = 0, along the path, shown as a dark (blue) 
line in Fig. 3. Eventually the Stark shift ceases, and the system point evolves in the vertical plane of As = 0. As the 
driving field diminishes, the system point moves towards the origin. Along this path there is an interaction O, and 
thus the adiabatic state $0 is a superposition of the two diabatic states [see eqn. (24)] . The population is thus shared 
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between the states tpi and ^3- We will show, in section VA, that the relative composition of this superposition can 
be altered by changing either the detuning A 2 or the ratio of the two Rabi frequencies Qj_ and 2 at the conclusion 
of the pulse sequence. 

An earlier paper [24] discussed similar adiabatic evolution, but with a delay between the two driving laser pulses. In 
that work the Stark-shifting pulse coincided with one of the driving fields [24] . Under these conditions, the system is 
divided between two adiabatic states of the Hamiltonian of eqn. (17). The relative phase of the resulting superposition 
depends on the time integral of the difference between the two associated eigenvalues. Many applications require the 
control of this dynamical phase, i.e. the time integral of the Rabi frequency must then be controlled. We have here 
considered the case where the driving pulses are synchronized and delayed with respect to the Stark pulse. In contrast 
to the sequence discussed in ref. [24], the dynamics of SACS follows a unique eigenstate and the relative phase of the 
superposition is reduced to the sum of the two driving frequencies. 

B. Optimization of adiabatic evolution 

As can be seen from viewing the energy surfaces of Fig. 3, many paths in the two-dimensional parameter space link 
the initial state with the desired final-state superposition. It is only required that the evolution be adiabatic and that 
the system remains on the surface associated with the adiabatic state $o- In general, the adiabaticity requirement 
implies that the time integrated Rabi frequencies are much larger than 1. This, together with the requirement that 
the transition to the intermediate state is near resonant, |A 2 | < fi max , leads to the condition 

A 2 T>1 (25) 

where T is the duration (FWHM) of the driving and Stark pulse, assuming Gaussian temporal shapes. As was done 
earlier [28], we fix the detuning A 2 and minimize the non-adiabatic losses by adjusting the peak values of the Rabi 
frequency fi max , the Stark-shift A™ ax and the delay between the pulses. That earlier work, modeling a two-level 
system driven by a chirped pulse, found that the non-adiabatic losses are minimized (for fixed pulse shape and peak 
Rabi frequency) when the system point follows a trajectory of constant energy (i.e. a level line) in the adiabatic 
energy space generated by the Rabi frequency and the detuning. 

In a temporal representation, the non-adiabatic couplings between two adiabatic eigenstates are minimized when 
the associated eigenvalues exhibit parallel temporal evolution. Following ref. [24], we used the criteria of parallel 
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evolution to find laser parameters that minimized the non-adiabatic couplings between the energy of the adiabatic 
state followed by the system and the closest alternative adiabatic energy. 




FIG. 4: (Colour online) Contour plot of the difference of the two lowest eigenenergies of the Hamiltonian (17) as a function of 
Q and As (in units of A2). The white line is an ideal trajectory that would optimize adiabaticity during population transfer. 
The dark trajectories are associated with the realisation of the transfer with sine-squared pulses. In contrast to the dotted 
(red) trajectory, the dashed (blue) trajectory follows approximately the optimal level line. 

Figure 4 shows the difference of the two lowest energy surfaces of Fig. 3. This contour plot exhibits level lines 
continuously connected to state ipi and to the desired superposition state. A path in the parameter space is represented 
by a closed loop starting and ending at fl = 0, A5 = 0. The white trajectory follows a level line. It is thus associated 
with a population transfer process which minimizes non-adiabatic losses. This trajectory will lead to parallel lines in 
the temporal evolution of eigenenergies for £1 ^ and A5 7^ 0. An inspection of the ideal loop gives the optimized 
peak amplitudes 

rr ax k 2.6A 2 , A£ ax k 1.5A 2 . (26) 

For a given temporal duration T, we first chose the value of A 2 satisfying (25). The conditions (26) give the values of 
the parameters r2 max T and A™ ax T which minimize the non-adiabatic losses. The delay between the pulses, and their 
shapes, have to be chosen such that the system follows the ideal trajectory as closely as possible. 

V. PROPOSED EXPERIMENTAL IMPLEMENTATION OF SACS 
A. SACS in mercury atoms 



We suggest a possible implementation of the proposed technique to the preparation of a coherent superposition that 
could enhance four-wave mixing in mercury vapor. The three states of interest are chosen from the degenerate energy 
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levels 6 1 S U , 6 3 Piand 7 1 5o, as indicated in Fig. 1. When the laser fields have definite polarizations, as we assume, the 
resulting magnetic-sublevel selection rule makes the excitation ladder a simple three-state system. Table I lists the 
relevant properties of these states. 



Transition 


A [nm] 


A12 [l^s- 1 ] 


d [10- 30 Cm] 


6 3 Pi -» 6^0 


253.728 


0.08 


2.15 


7 1 S -» 6 3 Pi 


407.898 


0.0088 


1.46 


6^1 -» 6^0 


184.950 


7.46 


12.94 


7 1 S - 6^1 


1014.254 


0.271 


31.67 


7 x Pi -» 6^0 


140.262 


0.226 


1.49 


7 x Pi -» 7 1 S 


1357.422 


0.150 


36.44 


9^1 -» 6^0 


125.056 


0.995 


2.63 


9^1 -» 7 1 S 


623.602 


0.017 


3.78 



TABLE I: Spectroscopic values of the considered transitions in mercury. A is the wavelength of the transition. Here A12 is the 
Einstein coefficient for the spontaneous emission and d is the dipole-transition moment. 

In the following simulations we assume laser pulses whose temporal shape follows a sine-squared pattern within 
one period. Such pulses are very similar to the Gaussian pulses often found in experiments, and they have the 
mathematical advantage of vanishing identically outside a finite pulse duration; they have finite support. The use 
of such pulses simplifies the identification of temporal regions where the two interactions vanish, and for which the 
adiabatic states have exact analytic expressions. Specifically, we used the profiles, 

fi(t) = O max sin ((* - t) IT) 2 for t G [t, ttT + t], = elsewhere, 

A s {t) = A^ ax sin(t/T) 2 for t e [0,nT], A s (f) = elsewhere, 

with equal widths of T = Ins. From table I we deduce the relation between the Rabi frequencies and the laser 
intensities as 

SIP [s _1 ] = 5.6 x 10 7 x ^/h [W/cm 2 ], f2™ ax [s^ 1 ] = 3.8 x 10 7 x [W/cm 2 ]. (27) 

For our simulations we chose the peak intensities as 1\ w 0.9 MW/cm 2 and I2 ~ 1.9 MW/cm 2 for the transitions 1 
and 2. These choices give equal peak Rabi frequencies of f2 max = 52.1 ns _1 for both transitions. The delay between 
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the driving pulses and the Stark pulse is r = 0.8 ns. The Stark shift of state 3 is evaluated from the expression 

As(t) -^ x y M m 

j J 

summed over all the non- resonant intermediate states j. Here £s and uus are the electric field envelope and the carrier 
frequency of the Stark pulse respectively. The energy of state 3 shifts by a maximum of A™ ax = 30.5 ns" 1 . This 
shift can be achieve using a Stark laser with the wavelength 1064 nm and with the intensity Is ~ 16 MW/cm 2 . We 
evaluate the Stark shift, taking into account the relevant states 6 1 Pi and 7 1 Pi, as 

As ax I 8 " 1 ] = -1885 x I s [W/cm 2 ]. (29) 

We neglect the shifts of states 1 and 2, as is usually the case for lower-lying states. This approximation was confirmed 
from estimations of the Stark shifts. 

These parameters have been chosen based on the following criteria. The pulse durations were fixed at T = 1 ns. 
The fields 1 and 2 were detuned by A 2 = — A 3 = 20/T to satisfy the condition (25). The intensities were chosen 
such that ft max = 2.6A 2 . The Stark shift satisfies A^ ax = 1.5A 2 to fulfill cqn. (26). The delay between the Stark 
and driving pulses was adjusted such that, in Fig. 4, the trajectory of the system point in the parameter space (blue 
dashed line) follows as closely as possible the ideal level line (white level line). At the beginning of the sequence, 
when the driving pulses 1 and 2 are still absent the populated adiabatic state coincides with the initial state Vi- At 
the end of the process, when the Stark pulse is absent, this adiabatic state corresponds with the dark state of eqn. 
(13). Because the Rabi frequencies f2i,f2 2 remain equal as they diminish, the populations Pi(t) and P^{t) are both 
equal to 0.50 and the state of the system, for t > tf, is 

(t) w -L U - e-^+^Vsl (30) 
V2 L J 

with uj\ + uj 2 = W12 = 1.20 x 10 16 s _1 , corresponding to a wavelength of 157 nm. 

The optimization of the pulse parameters to minimize the non-adiabatic losses can produce a variety of superposi- 
tions. By changing the relative intensity of the pulses, and hence the relative peak Rabi frequencies (fi max = 2.2A 2 , 
and A™ ax = 2.1A 2 ), we cause the trajectory of the system point to deviate from a contour line (red dotted line in 
Fig. 4) and the populations become Pi — 0.48 and P 2 = 0.52. 
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Figure 5 shows results of numerical simulation of SACS in mercury. The upper frames show the pulse sequence, 
while the middle frames show the time varying adiabatic eigenvalues. The thick line marks the eigenvalue associated 
with the adiabatic state that coincides with tpi before the pump pulses are present and with the dark adiabatic state 
$o, eq n - (13), after the Stark pulse vanishes. The lower frame shows the populations. In this example, they evolve 
into a 50:50 superposition of states 1 and 3. 




2 4 
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FIG. 5: Numerical simulation of SACS in mercury. Upper frame: Stark shift (dotted line) and driving (full line) Rabi frequency 
(fii = Q2 = fi) as a function of time. Middle frame: eigenvalues versus time. Lower frame: populations P„(t) as a function of 
time (Pi, dotted line; P2 full thin line and P3 full thick line) 

A maximum coherence, once formed, can be used in any subsequent frequency mixing process. The atoms act like 
a local oscillator at the two-photon transition frequency u)\ +ui 2 - This can beat with an additional radiation field of 
frequency W3 to generate new radiation fields, e.g. at the sum of all the contributing frequencies W4 = w\ + lo 2 + W3. 
When applied to four-wave mixing in mercury, the introduction of a tunable probe laser field with wavelength in the 
visible regime will generate radiation deep in the vacuum-ultraviolet spectrum. 

In contrast to conventional frequency mixing techniques, the maximum coherence maximizes the induced polariza- 
tion in the medium, and thereby enhances the efficiency of the frequency mixing process. Moreover, the coherence in 
the medium decays only as slowly as the excited state lifetime, and so the mixing process is also possible when the 
probe laser is delayed with respect to the driving fields. This feature has no counterpart in conventional nonlinear 
optics, which requires coincident radiation fields. 
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B. Alternative superpositions 

The proposed SACS technique offers the possibility of constructing any superposition of states 1 and 3 with con- 
trollable weights by changing the two-photon detuning A 2 + A 3 . Our analytical analysis and concomitant discussion 
assumed two-photon resonance. It is instructive to see the consequences of relaxing this condition. Figure 6 shows 
the population of states 1 and 3 after the interaction with the pulse sequence for various detunings A 2 + A 3 . The 
pulse parameters are identical to the ones of Fig. 5. When the two-photon detuning is far from zero and negative 
((A 2 + A 3 )T <C 0), the population is left in the ground state. In the opposite case, ((A 2 + A 3 )T 3> 0), there is 
complete transfer into state 3. For intermediate detunings the populations are shared between states 1 and 3 with 
various weigths. When two-photon detuning is present additional phases appear in the expression of the superposi- 
tion (see [27]). The relative weights of the components and ^3 can be altered by changing either the two-photon 
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V A 3 ( ns ~ 1 ) 

FIG. 6: Populations Pi (dotted line), and P3 (full line) at the end of the pulse sequence as a function of the two-photon 
detuning A2 + A3. 

detuning or the ratio of the Rabi frequencies f2i(i)/ri 2 (i) at the end of the pulse sequence. Indeed, the superposition 
state created by the two pulses is exactly the STIRAP dark state, eqn. (13). From the latter expression, it is obvious 
that the Rabi frequencies at the end of the interaction can be used as control parameters for the superposition of 
states tpi and ^3- 

C. Comparison with F-STIRAP 

The Stark pulse of SACS allows successful implementation with a large variety of pulse shapes. The dark state 
in F-STIRAP, eqn. (13), connects state ipi to the target superposition of eqn. (4), as long as the pump and Stokes 
pulses maintain a constant ratio of Rabi frequencies as they vanish. In SACS, the Stark pulse lifts the degeneracy 
of the adiabatic states, thereby permitting implementation with simultaneous driving pulses 1 and 2 of similar shape 
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and width. Figures 7 and 6 illustrates this flexibility with simulations using Gaussian pulses, 

n p (t) = r>™ ax e - ( ' /T)2 (31) 

Sl„(t) = n^e-^- T ^ T \ (32) 

A S (t) = A^x e -((t-r^l)/T) 2 _ (33) 

Each of these figures depicts the populations of states 1 and 3 at the end of the SACS process as a function of the 
Stokes peak amplitude and of the delay between the Stokes and pump pulses. For Fig. 7 the Stark shifting pulse 
is absent (the process thus corresponds to STIRAP). Population transfer exhibits very clear Rabi oscillations, as 
indicated by the curved bands. By contrast, the simulations of Fig. 6 include a Stark shift. The contour patterns 
here are indicative of adiabatic passage, as contrasted with Rabi oscillations. 




FIG. 7: (Color online) Contour plots of populations in state 1 (lower frame) and state 3 (upper frame) at the end of the pulse 
sequence as a function of the Stokes Rabi frequency (f2™ ax ) and of the delay (r) between pump and Stokes pulses for the 
F-STIRAP process. The Stark field is here absent. The parameters are the following: The detunings A2 = —A3 are equal to 
20/T. The pulse shapes are Gaussian, as specified in eqn. (33), with f2™ ax T = 40. The four outset plots show pulse sequences 
for four couples of values (f2™ ax ,r) correspond to the dots (a), (b), (c), and (d). In those plots the pump field Op is a full line 
and the Stokes field fls is a dashed line. 

Viewing Fig. 7 one can identify two distinct regimes for the F-STIRAP process. For a small delay between the 
pump and Stokes pulses, Rabi oscillations of the population occur between the states 1 and 3. For sufficiently large 
delay, the dynamics becomes adiabatic and there occurs an equal-weight superposition of states 1 and 3. 

Viewing Fig. 8 we see that when the Stark shift is sufficiently large the Rabi oscillations disappear and the dynamics 
becomes adiabatic even if there is no delay between pulses 1 and 2. The weights of the superposition components 
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FIG. 8: (Color online) Contour plots as in Fig. 7. Parameters are as in that figure, but there is a Stark shift, Ac? ax T = 50. In 
the outset frames the pump pulse Op is a full line, the Stokes pulse Qs is a dashed line and the Stark shift As is a dotted line. 

then depend on the ratio of Rabi frequencies Qi(t) / 0,2(t) at the end of the sequence. 

When a Stark shift of the energy of state 3 occurs the two-photon resonance condition is generally not fulfilled, 
and the adiabatic state followed by the statevector does not coincide with the dark state of STIRAP, eqn. (13). 
Under these circumstances the intermediate state ip2 is populated when both Stark and driving pulses are present. 
For adiabatic transfer to succeed it is necessary that the lifetime of the intermediate state be much longer than the 
duration of the overlapping between pulses. 



D. Comparison with two-photon half-SCRAP 

The SCRAP technique, applied to excite high-lying states, is usually based on a two-photon transition with no 
resonantly driven intermediate state. Thus it requires higher intensities than processes that rely on a combination 
of single-photon near-resonant excitations, as SACS does. To illustrate the main differences between SACS and half- 
SCRAP, we calculate the intensity of the driving field we would need to excite the transition 6 1 So-7 1 Sq by half-SCRAP 
with a similar Rabi frequency to the one used in the SACS technique. As above, the states 1 and 3 are taken to be 
6 1 Sq and 7 1 Sq. A pump laser at wavelength 313 nm drives a two-photon transition between these two states. The 
effective two-photon Rabi frequency is calculated as 

^=%XE j -E 1 J -tu p (34) 
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In this sum we included only the intermediate states 6 3 Pi, 6 1 P\, 7 1 P\ and ^Pi. Other states give only minor 
contributions. We calculate the effective two-photon Rabi frequency to be 

[s- 1 ] = 37 x I p [W/cm 2 ]. (35) 

The adiabatic condition S!£Jf x T m 50 3> 1 can be fulfilled for a pulse duration of T = 1 ns with a pump intensity 
of I p wl.3 GW/cm 2 . The pump intensity needed for half-SCRAP is thus three orders of magnitude larger than the 
intensities of the pulses 1 and 2 used in the SACS technique for the same adiabatic condition 0™ ax T = Vlf ayi T w 50. 
On the other hand, with the SCRAP technique the intermediate states are far from resonance, and so they are never 
populated during the dynamics. 

VI. CONCLUSIONS 

Four-wave mixing, and other nonlinear optics processes, are enhanced by preparing the nonlinear medium in a max- 
imally coherent quantum state, thereby maximizing the induced polarization field. For such purposes it is necessary 
to create superpositions of nondegenerate quantum states. We have introduced an alternative adiabatic technique, 
SACS, for the preparation of coherent superpositions of two nondegenerate states, 1 and 3, coupled by a two-photon 
process. The technique uses two pulses, each near resonance with an intermediate state 2. The time dependent 
Rabi frequencies of these two pulses are assumed to have identical pulse shapes. Additionally, a Stark-shifting pulse 
manipulates the energy of state 3. The weights of the components of the superposition can be controlled by adjusting 
the two-photon detuning or by changing the ratio of the driving pulses at the end of the sequence. We gave a set of 
optimal parameters that minimize the non-adiabatic losses. 

The SACS technique has potential advantages over alternative techniques for preparing such nondegenerate su- 
perpositions. The SACS method requires less laser power than half-SCRAP. It also offers a larger choice than does 
F-STIRAP for the pulse shape of the fields and for the delay between pulses. 

We simulated a possible implementation of the SACS process in mercury atoms. Our results suggest that it should 
be feasible to implement SACS and thereby to enhance four-wave mixing significantly. 
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